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\ Abstract: Large Rapidity Gap Events in Deep Inelastic Scattering are discussed in 

j^ I terms of lightcone wave functions for quarks and gluons inside the photon. It is shown that 



m 






this approach is consistent with earlier, conventional Feynman diagram calculations. An 
updated parametrization for the cross section is given and a numerical analysis presented. 



: 1 Introduction 

P^ ■ With the start of HERA a new tool has become available which allows a significant progress 

t^^ , in sorting out the longstanding puzzle about the nature of the Pomeron [|l| |^ . Theoretically 

^ ! the Pomeron is difficult to tackle. It was introduced phenomenologically as a simple moving 

1-^ I pole in the complex angular momentum plane (Pomeron trajectory) analogous to meson 

exchanges, however, compared to mesons there is no clear evidence for associated bound 

states in the s-channel (candidates are glueballs). The Pomeron intercept is slightly above 

^ . one which translates into a slowly growing total and elastic cross section in hadron reactions 

(soft Pomeron). Meson trajectories, on the other hand, are below one and give subleading 

K/i ; contributions at very high energies. Within the framework of perturbative QCD the Pomeron 

;h ' is associated with the resummation of leading logarithms in s (total energy) which results 

in a more complicated branch point singularity instead of a simple pole 0] (hard Pomeron). 

The major shortcoming of this leading log(s) approach is the ignorance of nonperturbative 

contributions and the neglect of unitarity corrections which are relevant for the complete 

formation of the Pomeron |]^, ^. Other approaches propose a combination of soft and hard 

Pomeron (see for example 0, ^) where the Pomeron intercept is controlled by the relevant 

scale of the process. 

Due to its zero color charge the Pomeron is associated with the occurrence of rapidity 
gaps at very high energies. The simplest form of rapidity gap events beside elastic scattering 
is single diffraction where only one of the incoming particles dissociates (the virtual photon 
in Deep Inelastic Scattering) whereas the other (the proton) stays intact. In most cases the 
proton escapes undetected, but it is surrounded by a Large Rapidity Gap, a fact which is 
used in experiment to define diffraction. In the following we will assume the proton not to 
decay. 
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Deep Inelastic Scattering exhibits the nice feature of having a small colorless particle, 
the virtual photon, in the intial state. We will make use of this fact and shift the focus 
from the less well defined Pomeron to the virtual photon which we believe is perturbatively 
calculable, i.e. its content of quarks and gluons can be determined. Going into the target 
rest frame the following picture emerges: a fast traveling photon dissociates far upstream the 
proton-target into a quark-antiquark pair which evolves into a more complex partonic system 
before the actual interaction takes place. The initial separation of the quark-antiquark pair 
in impact parameter space is of the order of 1/Q [Q"^ is the photon virtuality). The final 
partonic system, however, which takes part in the interaction covers a much larger area 
roughly the size of a hadron. It appears to be rather natural to refer to the leading order 
quark-antiquark pair as a color dipole, it is less trivial, though, for a multi particle state. 
Still, it may be shown that for leading twist contributions the concept of a color dipole can 
be extended to a more complicated final state. The argument goes as follows: hard QCD- 
radiation generates a bunch of partons strongly ordered in impact parameter space. The 
large distance is marked by the last quark or gluon in the chain of emissions whereas the 
remaining partons are confined in a small area (short distances). The separated parton on 
the one hand and the confined system of partons on the other hand form a new effective 
dipole. Both types, the quark dipole and the gluon dipole, can be described in terms of 
lightcone wave functions. The gluon dipole, although it is of higher order in perturbation 
theory, is of particular relevance when the invariant mass M of the diffractive final state 
becomes large. The dipole picture is certainly limited in its applicability, but it works fine 
when Q^ is the leading scale in the process, i.e. M'^/Q'^ is not extremely large. 

For the interaction of the color dipole with the proton we employ the two gluon model [0. 
All possible couplings have to be added up in order to retain gauge invariance. The two gluon 
model itself may not account for the full structure of the Pomeron, but it can be generalized 
to the exchange any number of gluons, since all gluons hooked on to one leg of the dipole can 
be merged into a single effective vertex. The strategy in this paper is to factorize the Pomeron 
structure from the dipole according to the /cj-factorization theorem [^. The dipole part is 
calculated whereas the Pomeron part acquires a phenomenological parametrization. A full 
QCD-treatment is not feasible at the present time. All free parameters will be determined 
from inclusive Deep Inelastic Scattering data {F2) and then used for diffraction. 

In the following section we will introduce the lightcone wave function formalism which 
includes an improved expression for the gluon dipole. In section 3 the model for the Pomeron 
is specified and a fit to F2-data which determines all parameters is performed. The Pomeron 
model is then combined with the lightcone wave functions and the cross section for diffraction 
in Deep Inelastic Scattering is calculated. It is shown that these results are consistent with 
earlier approaches based on more conventional calculations using Feynman diagrams (section 
4). In section 5 the cross section is numerically evaluated and the main results are discussed. 



2 Light cone Wave Functions 

One should take the terminology 'wave function' not too literally, since the state it is meant 
to represent is not real but virtual. Still, it provides an effective and intuitive description 
of quark dipoles and gluon dipoles inside the photon in close analogy to a quark-antiquark 
boundstate (quarkonium). An unpleasant property which can be traced back to the pointlike 
nature of the photon is the lack of normalizability of the wave function unlike for real states. 
In order to find the correct normalization one has to go back to the corresponding Feynman 
diagrams 
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The nice feature is that once having determined the lightcone wave 
function one can easily study single gluon and multi gluon exchange, although the wave 
function formalism for the gluon dipole is restricted to color zero exchange, i.e. it works for 
two or multi gluon exchange. 



q — k 
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Figure 1: Quark dipole (a) and gluon dipole (b) 

The virtual photon can be transverse polarized (transverse with respect to the lightcone 
vectors q' = q + XbP and p, q is the photon momentum and p the momentum of the incoming 
proton) as well as longitudinally polarized. The two helicity states 7 = ±1 follow from the 
projection on the transverse vectors (l,i) and (1,-i). For the left and right handed quarks 
we introduce the quark helicity h = ±1. The quark momentum k may be parametrized like 
k = aq'+(3kP+kt (Sudakov parametrization) which simultaneously fixes the decomposition of 
the antiquark momentum: q — k = {l—a)q'+{—XB — P)p — kt. Using complex notation for the 
two dimensional vector kt we may write the lightcone wave function for the quark-antiquark 



state as: 
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where Q^ = —q^ is the virtuahty of the photon. Eq.(^ shows the hghtcone wave function of 
the longitudinally polarized photon. 

Before discussing the properties of the lightcone wave function it is necessary to have a 
closer look at the kinematics. We assume that either the quark is offshell and the antiquark 
onshell or vice versa (here: (g — fc)^ = and k'^ ^ 0, fig.|]). In the frame that we choose the 
photon moves fast and the two quarks roughly carry the momenta ag' and (1 — a)g' whereas 
the other components are small. Any subsequent high energy scattering does not change the 
a-component. The offshell quark with the momentum k becomes onshell after receiving a 
small fraction xp of momentum along p while being scattered. The momentum of the quark 
changes from k to k with k^ = 0. Using the mass-shell condition one finds: 
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W is the total hadronic mass {W'^ = 2q' ■ p). The missing mass M is simply given as the 
total energy of the two outgoing quarks: 
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We now include the emission of a gluon in our discussion. At large M the gluon is 
well separated in rapidity from the gg-pair and becomes the dominant configuration over 
the exclusive gg-pair production. The latter is suppressed by a power in M^ (spin 1/2- 
exchange). Unfortunately, the three particle Fock state is much more complicated, and a 
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rigorous construction of the wave function which is consistent with Feynman rules and vahd 
for all kinematics has not been achieved, yet. A simplification occurs when only the leading 
twist and leading log(Q^) contribution is considered. In this case the distance in the impact 
parameter space between the quark and the antiquark is much smaller than the distance 
between the quarks and the gluon. The gg-pair on the one side and the gluon on the other 
side form an effective color dipole similar to the exclusive gg-pair that we have considered 
before (fig.|l|.b). 

For the gluon dipole we find the following wave function (after introducing the vector 
notation kt = [k], kf)): 



In view of the two-vector-particles state (an effective two gluon state) it appears rather 
natural to find a tensor representation for the wave function. In the triple Regge limit 
(TRL) with M^ much larger than Q^ the term a{l — a)Q'^ in the denominator of expr.(|^) 
may be neglected and only kf remains: 

^^TnM^k,) = , ' _^ '-^ . (6) 



a(l-a)Q2 k, 



The S-teiia which at first sight should be kept was also removed, since it does not depend 
on kt and drops out in any application due to subtractions (see eq.(p!9[) of section 4). The 
simple structure of eq.(^ was found earlier in refs. [[TI|, |13|. Expr.(|^), on the other hand, is 



valid for all masses and provides the natural extension of the Triple Regge result. 

We have to point out that the wave function introduced in eq.(|^) does not reproduce the 
amplitude for the single gluon t-channel exchange. Especially those contributions which are 
singular in the limit M ^ are absent. The reason for that lies in the fact that in the color 
singlet configuration all singular, soft terms cancel out. This cancellation can be illustrated 
by considering final state radiation off the gg-pair. The soft terms add up when the gg-pair 
is colored, they cancel each other, however, when the state is colorless. 

We have written eq.(^ in a symmetric way with respect to a and 1 — a in order to 
stress the similarity with the gg-dipole. The gluon dipole considered here is actually very 
asymmetric which is a consequence of the leading twist approximation where the internal 
virtualities are much smaller than Q^. From eq.(^ we conclude that a has to be much 
smaller than 1 in order to fulfill the condition k"^ <^ Q^. So, one could have set 1 — a in 
eq.(H) equal to 1 without reducing the accuracy of the formula. 

At very large masses M or small /3 the dipole picture becomes insufficient. Instead of 
logarithms in Q^ we have to sum up logs in 1/(3 or M"^ /Q"^. This leads to a new four-gluon 
t-channel state P, Q and a new singularity 1 + ti;4 in the complex angular momentum plane 



((l//5)'^*). But in contrast to conventional QCD-scaling violation the strong rise with 1//3 
does not imply a rise with Q^. 



3 Modeling the Structure Function F2 

We compute the structure function F2 at very low Xb when single gluon exchange gives the 
leading contribution. In contrast to the usual approach (gluon-boson fusion) the gluon is 
not onshell and needs to be described by a distribution over both, longitudinal and trans- 
verse, phase space components (fc^- factorization). Such a factorized distribution (JF) con- 
tains perturbative as well as nonperturbative contributions, and the aim is to find a suitable 
parametrization which gives a reasonable description of all low xb and low Q^ data. 




Figure 2: Inclusive Deep Inelastic Scattering (-F2). 

The fct-factorization theorem is the high energy or small Xb counterpart of the conven- 
tional (collinear) factorization theorem. The latter, when calculating the structure function 
F2, requires a convolution of the gluon structure function and the quark box with respect to 
the longitudinal momentum fraction whereas in the small xb regime the convolution with 
respect to the transverse momentum (kt) is more appropriate. At zero momentum transfer 
the Pomeron (in Deep Inelastic Scattering) is essentially the same as the unintegrated gluon 
structure function, and by fitting the F2-data one determines the Pomeron intercept ajp. In 
Deep Inelastic Scattering, however, the Pomeron intercept varies (depending on Q"^) rather 
than being a fixed number as in soft processes {ajp = 1.085). 

As we already mentioned one of the virtues of the wave function formalism is that we 
can use it for single gluon exchange (transverse momentum It) as well (see fig.^]). We find 



for F2 the following expression: 
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In the limit of large Q^ (at fixed but small Xb) and taking only the leading log(Q^) contri- 
bution the gluon distribution factorizes and the Zf-integral can be taken: 
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g{xB,Q'^) represents the conventional gluon distribution and J-'/ag is usually referred to as 
unintegrated gluon structure function. Introducing the Feynman parameter x eq.(|^) reduces 
to 
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This representation (see also |T6|, |17|) serves as starting point for further numerical evaluation. 
As ansatz for T we choose: 
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Qo is set to IGeV (proton mass), and since only small xb are considered we introduced 
xq = 0.05 as normalization point. The Pomeron intercept has the following parametrization: 



B ln[ln(QVg2) + 1] if Q^ > q2 

if g' < Ql 



ajp{Q') = 0.085 + <i ;; "^^"^^^ '^0^ ' ^^ :; ;^2<5°2 (n) 

A soft Pomeron {ajp = 1.085) intercept is assumed when the scale Q^ falls below Q^. This 



behavior has experimental support form the BPC-data |T8 



The ansatz in eq. (|lO|) has the following two basic ingredients: first, JF scales like 1/Qq for 
It = 0, i.e. an effective cutoff at the scale of Qq is introduced which eliminates the singularity 
related to the gluon propagator l/l^. The scale Qo roughly represents the inverse size of the 
hadron and because hadrons are colorless all gluons with a wave length larger than the 
size of the hadron decouple. One important consequence is the vanishing of the structure 
function when Q^ approaches zero. The second ingredient is a scale (Q^-) dependent Pomeron 



intercept. It takes care of the experimental fact that the small-xs rise becomes weaker when 
the scale decreases and and finally turns into the soft behavior below Qq. The parameter 
B in (0) has to be determined from data. Two more parameters come along with B (see 
eq.(|lOD), A fixes the absolute normalization and C corrects the strong scaling violation which 
results from the scale dependent parametrization of the Pomeron intercept. 



The data for the fit are taken from HERA |18, 19| below Q^ = 50GeV^ and from E665 |20 



including the smallest Q^-values. The fit gives the following values for the three parameters 
A, B and C: 
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These parameters will be used in the following to predict the diffractive structure function 



4 Diffraction 

For the diffractive cross section we use the same conventions as for the inclusive cross section 





a) 



b) 



Figure 3: Wave functions with shifted argument. 

in Deep Inelastic Scattering, i.e. we decompose it into a transverse and a longitudinal part 
according to the different polarizations of the virtual photon: 

da 
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(13) 



where Wt and Wi are the transverse and longitudinal projection of the hadronic tensor W^'^: 
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with s = 2q' ■ p and q' = q + XbP- The transverse tensor g^'^ is defined as 
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(14) 



(15) 



The momentum transfer is set to zero, a restriction, which is justified by the fact that the 
cross section peaks at t = 0. 

At leading order we have to consider the coupling of two t-channel gluons to a quark- 
antiquark pair. The four possible couplings are represented by four terms involving the 
wave function with shifted transverse momenta (see fig.^a). These shifts correspond to the 
transverse momenta carried by the gluons {It and —It) when the total momentum transfer is 
zero: 
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The last expression results after integration over the azimuthal angle. 

For the total contribution to Wt and Wi one has to take the square of the amplitude, 
i.e. basically the square of the previous expression with summation over the corresponding 
helicities: 
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We have here substituted k^ using eq. 



in combination with (3 = Q'^/{M'^ + Q^). Similar 



expressions can be found in [^, |TT|, [12 



As in the previous section we have absorbed the strong 
couphng constant into J-". This way the remaining expressions become free of parameters. 

By taking the hmits /5 — >■ 1 and /3 — *> one can study the main properties of eqs.(|l^) 
and (plsl) . For /3 — >• 1 eq.([T7|) vanishes proportional to (1 — /?) whereas eq.(|l^) gives a 
finite contribution. Hence, at small masses the longitudinal contribution is larger than the 
transverse contribution. A similar observation is made for vector meson production in Deep 
Inelastic Scattering. One, however, has to be aware that the longitudinal part is of higher 
twist, so that the ratio of the transverse part to longitudinal does not vanish like M^/Q^ at 
large Q^, but is roughly of the order of M'^/Q'^ with a logarithmic enhancement (see also 
ref.[^). Taking the second limit, /5 — > 0, one observes that the transverse contribution is 
finite, i.e. in terms of the mass M: the cross section vanishes like 1/M^ as expected for a 
spin 1/2-exchange. The longitudinal part ( jTSf ) has an asymptotic behaviour proportional to 
/3^, i.e. is negligible at large masses. 

For the configuration with a gluon in the final state we go back to eq.@. As in eq.(|l^) 
we encounter four terms (see fig.^.b): 
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To get the complete contribution for Wt one has to, first, take the square of the previous 
expression and, second, add the contribution for the perturbative splitting of a gluon into 
two quarks: 
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Two new variables were introduced in this equation, the longitudinal momentum fraction z 
(relative to the Pomeron momentum) and the virtuality /c^ of the t-channel gluon which is 
connected to the quark-box. }^ is related to the transverse momentum by fc^ = fc^/(l — 2;). 
The previously used variable a was substituted by means of the equation a(l — a)Q^ = z)^ 
where we assume that a -C 1 according to the leading log(Q^) approach adopted here. A 
factor of two arises in the integration from the opposite and symmetric limit (1 — a) ^ 1. A 
contribution to Wi is negligible in this case. We, again, take the two limits /5 — > and /3 — > 1 
in order to understand the basic behavior. For small /3 the region of small z dominates, 
so that the second line in eq.(^OD can be approximated by setting z to zero. Integrating 
over z then results in 1//?, i.e. the cross section is divergent when /3 approaches zero. The 
small /3 or triple Regge hmit has already been consider before in refs. [|], 0, |l^, [13|, |15| , and 



the results are found to be consistent with our calculations. One can also start with '^^^i^ 
(eq.(|])) which directly yields the triple Regge limit of eq . (|20|) . Taking the opposite limit 
/? ^ 1 which was not covered by previous calculations one finds that eq.(PD|) vanishes like 



1^1 — /?)^. This result emerges by taking the limit z ^ 1 resulting in {1 — z)"^ and integrating 



over z. A derivation of eq.(^) based on Feynman diagrams is given in ref. [ 

It is important to note that the function JF is assumed to be universal and should be the 

same for all three eqs.(|T7|), ([T8| ) and (|20|). The only parameter that enters in eq. (|20|) is the 

strong coupling constant a^. It depends on the scale somewhere between Qo and Q. 

One may introduce a simplification by taking the limit /< — > in the wave functions of 

eqs.(0) and (|19|). This limit gives the leading contribution provided that kf/{l — /?) is much 

larger than Ql: 
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One may call this result leading log{k'^/Ql) approach for the remaining integration over If 
(eq.(|])) is logarithmic (one propagator 1//^ is hidden in JF). It also indicates the limitation 
to large k'^. 
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We again convert the previous results into contributions to Wt and Wi. Substituting a 
by means of eq.(|^) and introducing the virtuahty k"^ = kf/{l — (3) we find: 

- EQ?f 4/3^(1-/3) /J'^ [a^xjpgix^^k^Ql)}' (22) 
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Eq.(|2^) was derived earher in refs. [p^, p^ |25| , p6| and eq.(^) in ref. |2^. The virtue of 



this approach is the fact that we do not need the unintegrated structure function, instead 
we can use the conventional structure function. The shortcomings, however, are first of all 
the need of a cutoff k"^ which can only be realized by requiring jets in the final state. Second, 
one recognizes that in ( ^2]) the small (3 region is strongly suppressed and not constant as 
anticipated earlier. In this regime next-to-leading log(fc^/(5g) corrections become important. 
Corrections of this type have been explicitly calculated in ||2^. We face a similar situation 
for the longitudinal contribution eq. (|23|) which is zero at /5 = 1/2. Again next-to- leading log 
contributions become relevant here. 

To complete the discussion on the leading log(/c^/(5o) approach we give the corresponding 
formula for gluon production: 
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As before we take the square of the previous expression and rewrite the result in terms of 
Wt . The procedure is similar to the derivation of eq. (|20|) where a{l — a)Q'^ is substituted 
by zk"^ and the splitting function for gluons into quarks is added: 
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This result has been derived earlier in refs. [p^, yA > |25| by direct calculation of Feynman 
diagrams. 



5 Numerical Results 

With the formulae for Diffraction at hand and JF determined by the inclusive F2 (see section 
2) we can numerically evaluate F2 ■ To this end we note that F2 and F^ are related to Wt 
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and Wi by 
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^ xjp V 47r 27r ' ^ ' 

' xip Air 



{[l + il-yy]F,''{P,Q^xjp)-2x,y'Fi''{f3,Q\xjp)} . (27) 



The longitudinal structure function Fi is usually negligible due to the accompanying factor 
y'^ which is experimentally small in most cases (y is the energy loss of the electron). But the 
longitudinal contribution is not completely lost, since it still appears in F^ (see eq. (p6D ) 
We insert the formfactor l/(/^ + Ql) which belongs to JF into eqs.(0), (|TB|) and 



and perform the integration over l^ analytically. For convenience we introduce the variable 
V which is defined as v = Ql/^k"^ + q;(1 — a)Q'^) = Ql{l — P)/k^ for the quark dipole and 
■^ = Qo/^^ — Qo(l ~" ^)/f^t fo^ the gluon dipole. F^ is then presented in three separate 
contributions {Ff = Fa + Fb + FJ according to eqs.(^)(Fj, (pD(Fb) and (13) (FJ: 



Fa{l3,Q',x,) = —1-^ r^ /^ ^ G\x„lM^il-2^) . (28) 



:i-2/5)ln(i^ + I , ^ ^^^"^ -1 + 2/^1 ln(t;) 



61-/5 [ \I3) V^t;2 + 2(l-2/5)t; + l 



l-2(3 + v , Jv^ + 2{l-2P)v + l -{l~2P + v) 

In 



^v"^ + 2(1 - 2f3)v + 1 \y^t;2 + 2(l-2/3)?;+l + (1 - 2/3)i; + 1^ 
F,(/?,Q^x,) = — i— r^_^^i=G^(x„l/t;)^ • (29) 



QBoQ"^ JiPQl/Q^ ^3 /^ _ ^QlP X^ 



•H'°wn /.-.2(i-2,>.r ^|'°'"' 



, >2 _^ 2(1 - 2/3)t; + 1 -{1-2/3 + v) 
In 



t;2 + 2(1 - 2P)v + 1 V \P^ + 2(1 - 2/5)t; + 1 + (1 - 2p)v + 1^ 



and 
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— m 




4 1 



/3^ 



'l+v-2z{l-z)] In - 




ln(t>) 



+ 




Since the measurement is not performed at t = we have assumed a simple exponential 
behavior in t, exp(i?Bt), with the slope parameter Bd taken from experiment [^ (fi^i = 
S.O/Gel^^). The integration over t leads to the extra factor l/B^,. For a^ we estimate a 
value of 0.25 which is a reasonable estimation for scales around 2 — 3GeV. The function 
G is defined in eq.([TO|) and enters the equations above without changing the parameters. 
In fig.| (first plot) we show the /^-distribution for fixed xp = 5.0 • 10""^ and Q^ = lOGeV^ 
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Figure 4: /3-spectrum. 

with separate curves for each of the three contributions F^, F\, and Fc. As was already 
argued analytically we find three distinct regimes in the /3-spectrum: i) small /3 where the 
configuration with a gluon in the final state dominates (-Fc), ii) medium (5 where the exclusive 
quark- ant iquark production with transverse polarization is dominant [Fa] and iii) large /? 
where the longitudinal production of quark-antiquark pairs takes over (-Fb). The second 



14 



10"' 



10" 



10'' 




(3=0.5 




Q = 1 GeV' 



10" _ 



1 O'" 1 0"' 



^ = 0.01 
(3 = 0.5 
(3 = 0.9 



Q =10GeV' 



(3=0.01 - 
(3=0.5 
(3=0.9 - 



Q =100GeV' 



1 0'' 1 0"' 1 0' 



|3 = 0.01 
(3 = 0.5 
/3 = 0.9 



lO"' 10"' 10"" 10" 



Figure 5: Xjp-distribution. 

plot displays the change in the shape of the /3-distribution with increasing Q^. It is rather 
flat around Q^ = lOGeV^ before it starts tilting when Q^ is further increased (higher twist 
suppression at the large end and logarithmic enhancement at the low end of (3). 

The next figure (fig.§) shows the xjp-distribution for fixed /3 and Q^ (values as indicated 
in each graph). The first row of plots starts at low j3 with rather fiat distributions in xjp 
(the slope is 0.17 due to soft contributions, ajp = 1.085). The distributions become slightly 
more curved and steeper when we move to larger /3. The change of the shape with /3 results 
from the nonfactorizing ansatz where the Pomeron intercept depends on an intermediated 
(variable) scale related to the size of the dipole {Ql/v in eqs.(p8D, ( pQ]) and (pOD). At low 
(3 this scale is small (large dipole) and the Pomeron is dominantly soft. At high j3 and in 
particular when the higher twist (longitudinal) part takes over the intermediate scale is, in 
average, rather hard (approximately Q^/4) which leads to steeper distributions. This effect 
can roughly be interpreted as large mass states having a larger radius than small mass states. 
The curvature in the double log plots is due to smearing when the scale is integrated. In the 
second row we see a similar behavior by changing Q^. At small Q^ the intercept is frozen at 
a low (soft) value close to 1.085. There is only a little, barely visible effect due to smearing. 
When Q^ is increased the intermediate scale is pulled up and one finds again a steeper and 
slightly curved distribution. In total we note that Regge-type factorization is violated, i.e. 
the Pomeron intercept depends on f3 and Q^. 

The scaling properties in diffraction are of special interest because they provide direct 
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Figure 6: Q^-distribution. 

information about the Pomeron structure. Fig.|^ shows the dependence of F:[' on Q^ for three 
different values of f3 (0.01, 0. 5 and 0.9) which as we know from fig.^ also distinguish between 
the three contributions Fc, Fa and Fb. (For very low Q^ (close to IGeV^) our formula for Fc 
has to be taken with care, since it is only computed to leading log{Q^) accuracy.) The value 
for xip is again fixed at 5.0 ■ 10~^. 

The prediction for the slope in Q^ is according to usual Q^-evolution negative at large 
/?. Fig.p, however, shows a rise at low Q^ for any value of (3, even for the longitudinal 
contribution {(3 = 0.9). The latter develops a maximum around Q^ = lOGeV^ before the 
asymptotic regime is reached and the l/Q"^ (higher twist) suppression sets in. For /3 = 0.5 we 
also see an increase with Q^ which then flattens out towards a constant (leading twist scaling) 
behavior. A rise over a certain range in the Q^-distribution is not completely surprising, 
since F:^ vanishes when Q^ approaches zero. What is surprising is the delay with which the 
asymptotic (scaling) regime sets in. This effect seems to be model dependent. Taking the 
hard Pomeron approach from ref. [^ as example (plots can be found in ref.[^]) the delay 
is even more pronounced. A precise measurement of the Q^-scaling behavior seems to be a 
promising tool to discriminate various Pomeron models. 

A Q^-scaling violation for rather large (3 which persists far into the asymptotic region 
can presumably not be reconciled with the dipole approach. An alternative scenario based 
on the hard component of the soft Pomeron would predict a log((5^)-type behavior [PT|. 



We have not presented a comparison with data here. This can, however, be found in [p2 



The theoretical curve in |3^ is based on the same model as presented in this paper, only the 
values for the parameters have changed slightly with little impact on the xjp-spectra. 
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6 Summary 

We have derived the cross section for diffraction in Deep Inelastic Scattering starting from 
two types of hght cone wave functions, one for a quark dipole (eq.(|^)) and the second for 
a gluon dipole (eq.(§)). The latter is of higher order in perturbation theory, since a direct 
coupling of photons to gluons is lacking. We have shown how the color dipole approach 
works for a multi parton state provided one stays within the realm of leading twist and 
leading log(Q^) accuracy (strong ordering in impact parameter space). Only the leading 
order quark- ant iquark pair forms a dipole for any kinematics. The light cone wave function 
formalism was proven to be consistent with Feynman diagram calculations, but it should not 
be confused with the general dipole approach of ref.|3^. As soon as the strong ordering in 



impact parameter is lost multiple dipoles may occur. The wave function formalism presented 
here is not able to cope with this configuration. 

The expression for the quark dipole (||) is rather well known, the second expression (^) 
for the two gluon dipole is new. It is valid over the complete range of invariant mass of 
the two gluons and therefore an extension of an earlier derived version which was limited to 
large masses (triple Regge limit, eq.(|D) 

As model for the Pomeron we have considered color zero two gluon exchange which is 
easily generalized to multigluon exchange. All gluons at each leg of the dipole merge into a 
single vertex, i.e. there is only a single interaction point in impact parameter space. We can 
factorize the dipole from the target (fc^-factorization scheme) and parametrize all unknown 
contribution in terms of an unintegrated gluon structure function. This factorization works 



for two simple perturbative gluons, for shadowing corrections ^^ and even for scattering in 



a nonperturbative classical field. It should be possible to reformulate the results in ref. p5 



along the line of our dipole approach. In the semiclassical approach of ref . [^ the gluon den- 
sity is directly related to the unintegrated structure function, and the Landshoff-Nachtmann 
model |]3^ can as well be identified with an appropriate unintegrated gluon structure func- 



tion. 

In this paper, however, a more phenomenological ansatz was chosen. A parametrization 
for the unintegrated structure function was established with parameters determined from 
inclusive Deep Inelastic Scattering {F2) and then inserted into the corresponding expression 
for diffraction. An important feature of our parametrization is the scale-dependence of the 
Pomeron intercept which results in a variation of the xjp-slope. The scale is roughly the 
inverse size of the dipole (not Q"^) which has to be integrated over. Its average is close to 
the soft scale of IGeV, but increases slightly with /3 and Q^ which then causes the xjp- 
distribution to become steeper. The lower limit for the slope is given by the soft Pomeron 
intercept. The /3-spectrum is subdivided into three regions each being dominated by the 
following contributions: i) the gluon dipole at small (3, ii) the quark dipole with transverse 
polarized photons at medium /3 and iii) the quark dipole with longitudinal polarized photons 
(higher twist) at large (3. Around Q^ = lOGeV^ the total spectrum is rather flat, but it starts 
tilting when Q^ is increased (falling from /? = to /? = 1). 

The Q^-distribution is of special interest because it helps revealing the structure of the 
Pomeron. We find that for jS > 0.5 the slope in Q^ is positive up to Q^ ~ lOGeV^ and then 
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flattens out (leading twist) or turns down (higher twist at large f3) (see fig.0). A comparison 
with data from HI was performed in ref.[^ where the agreement is found to be reasonable. 
The main deviation between theory and data is due to secondary exchanges which have not 
been included in this paper. The LPS-data from ZEUS p8| also seem to agree quite well. 

So far only leading order and most important next-to-leading order contributions have 
been taken into account. To obtain more precise prediction for j3- and Q^-distributions one 
needs to perform a complete Q^-evolution which will be subject of another publication. Also 
of interest is a next-to-leading order diffractive jet-analysis. This requires, however, a full and 
consistent next-to-leading order calculation which goes beyond the light cone wave function 
approach of this paper. 
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